Domain wall dynamics in soft magnetic materials 
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The dynamics of a domain wall in magnetostrictive materials is investigated. The domain wall 
is modeled by a d-dimensional interface moving in a d + 1-dimensional environment. Long-range 
demagnetization effects and quenched disorder are considered, while the external magnetic field is 
increased at constant rate. Exact expressions for the average interface velocity and susceptibility are 
obtained, resulting that the system is critical when the demagnetization constant and the average 
interface velocity vanishes. The critical exponents are computed to 0(e = 4— d). Our predictions are 
compared with numerical simulations and Barkhausen noise measurements reported in the literature. 
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Domain walls in fcrromagnets move following an irreg- 
ular motion in response to changes in an applied external 
magnetic field, leading to discrete jumps in the magne- 
tization. This phenomenon, known as the Barkhausen 
effect, has received a renewed theoretical and exper- 
imental interest, because of its connection with non- 
equilibrium critical phenomena like interface dcpinning 
and self-organized criticality 

Urbach et al Q have analyzed the role of demagneti- 
zation effects in the avalanche statistics. From numerical 
simulations and Barkhausen noise measurements, they 
conclude that the addition of an " infinite-range" demag- 
netization field results in self-organized criticality, while a 
local demagnetization field does not. On the other hand, 
Narayan [|2| has shown that long-range dipolar interac- 
tions give rise to self-similar dynamics. More recently, 
Cizeau et al || derived an equation of motion for the dy- 
namics of a ferromagnetic domain wall, including both 
long-range demagnetization effects and dipolar interac- 
tions. From numerical simulations, they conclude that 
long-range dipolar interactions decreases the upper crit- 
ical dimension from d uc — 4 to d uc = 2 ||. Moreover, 
they have also demonstrated that the demagnetization 
factor acts as a control parameter, criticality is obtained 
after it is fine tunned to zero. 

In magnetostrictive materials magneto-elastic inter- 
actions enhance the surface tension effect over dipolar 
interactions and, therefore, scaling exponents for two- 
dimensional domain walls should be lower than the mean- 
field ones. Experimental measurements of Barkhausen 
noise with different magnetostrictive materials corrobo- 
rate this conclusion E|,[|,|]l. 

In a recent work B we have analyzed the dynamics of 



a d-dimensional interface under a driving force increas- 
ing linearly with time, a restored force proportional to 
the local interface height H(x,t), and quenched disor- 
der (This model is equivalent to the Urbach et al model 
with local demagnetization effects). Extending previous 
renormalization group (RG) calculations developed for 
the constant force case [||J^] we obtained the dynamic z 
and roughness £ exponents in first order of e = 4— d. Here 
we use a similar approach to analyze the same model but 
with long-range demagnetization effects. The starting 
equation of motion is a particular case of that obtained 
by Cizeau et al || , excluding long-range dipolar interac- 
tions. It is demonstrated that the model with long-range 
demagnetization effects belongs to the same universal- 
ity class as that of local demagnetization effects . Our 
predictions are compared with numerical simulations and 
Barkhausen noise measurements, reported in the litera- 
ture. 

In the central part of the hysteresis loop, around the 
coercitive field, the magnetization process is mainly due 
to domain wall motion, while domain nucleation and ro- 
tation are absent. Moreover in most ferromagnetic mate- 
rials, due to magnetocrystalline forces or sample shape, 
the magnetization has preferred directions. In the sim- 
plest situation, there is a single easy axis of magnetiza- 
tion and the domains are separated by surfaces, parallel 
to the magnetization and spanning the sample from end 
to end. In this situation the domain wall dynamics may 
be modeled as the motion of a d-dimensional interface 
in a d + 1 environment, with height H(x,t). In magne- 
tostrictive materials the long range dipolar interactions 
can be neglected, resulting the equation of motion 



TV 2 H(x, t)+ht-e 



d d x' 



H(x',t)+T)[x,H(x, t)}, 
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where T is the surface tension of the wall, A represents 
a friction coefficient, ht is the applied field which is sup- 
posed to increase at rate h, L is the linear size of the sys- 
tem, and e is the demagnetization factor. Here rj(x, H) 
is a Gaussian uncorrelated noise due to impurities, with 
zero mean and 



{q(x, H)r)(x', H')) = S d (x - x')A{H - H'), 



(2) 



where A (if) is a monotonically decreasing even function, 
with a fast decay to zero beyond a distance a± . In eq. (|l|) 
we have explicitly taken into account that the demagne- 
tization constant scales as the inverse of the system size, 
which for an hyper-cubic lattice scales as L d . Hence, the 



demagnetization factor e does not depend on system size. 
Eq. (0) is similar to Edwards- Wilkinson equation with 
quenched noise, but the external field increases at rate h 
and there is an additional term due to demagnetization 
effects. This difference carry as a consequence that the 
interface is never pinned by impurities, but always moves 
with a finite average velocity v. 

A perturbative solution of eq. (|l|) can thus be found 
expanding H(x, t) around the flat co-moving interface vt. 
Taking H(x, t) = vt + w(x, t), with (w(x, t)) — 0, we ob- 
tain the following equation for w(x, t) 



TV 2 w(x , t) + (h- ev)t - \v 



d d x' _ 

——j-w(x, t) + r)[x, vt + w(x, t)}. 



(3) 



The average velocity is obtained using the constraint 
(w(x,t)) — 0, while fluctuations around the average 
will be analyzed using RG transformations. For these 
purposes is better to work with the equation for the 
Fourier transform of w(x,t), w(k,u>) 6J. The external 
field (h — ev)t gives a singular term of the order of uj~ 2 . 
This singular term predominates in the low frequency 
limit resulting, after imposing (w(k, u>)) = 0, 



h 

v = — . 

e 



(4) 



Another exact result can be obtained if one com- 
putes the low-frequency and long-wavelength suscepti- 
bility. Adding a source term tp(k, to) to the equation for 
w(k,oj) one obtains the generalized response function |j] 



1 



where 



[G (k, uj)}- 1 = Tk 2 + iXuj + eS(k, L). 



(5) 



(6) 



is the bare correlator and E(fe,w) is the "self-energy". 
8(k,L) is the Fourier transform of L~ d . In the thermo- 
dynamic limit (L — > oo) S(k, L) w I for k — > and zero 
otherwise. Since S(0, 0) = and Go(0, 0) = e^ 1 it results 
that the low-frequency and long-wavelength susceptibil- 
ity (or simply the susceptibility) is given by 



X = G(0,0) 



(7) 



This result is also exact. 

In the limit e -* the susceptibility diverges as x ~ 
e~ 7 with 

7=1. (8) 



The system is thus critical for e — > 0, i.e. when demag- 
netization effects becomes negligible. However, different 
behaviors can be observed depending on the value of h. 
If h is finite and we take e — > then, according to eq. (Q), 
the domain wall velocity will become infinitely large and, 
therefore, the noise term in eq. (||) will be reduced to an 
annealed noise. Thus, for (h > 0, e — > 0) the interface 
dynamics is described by the EW equation with annealed 
noise, which has an upper critical dimension 2. Hence, 
for d > 2 the interface will moves with a velocity given by 
eq. (||) and fluctuations around the average can be ne- 
glected. This picture corresponds with a suppercritical 
regime, where the external magnetic field predominates 
over disorder. On the contrary, for (h — > 0, e > 0) the 
system will be in a subcritical regime. The domain wall 
will moves with an infinitely small velocity and, there- 
fore, the quenched noise will lead to fluctuations around 
the flat co-moving interface, but these fluctuations will 
be not critical due to demagnetization effects. This fact 
becomes evident with the existence of a finite susceptibil- 
ity. Finally, the critical state corresponds with the double 
limit e — > and v = h/e — ► 0. 

To go further we perform a RG analysis of the problem. 
We integrate out the degrees of freedom in a momentum 
shell near the cutoff A and rescale k — > 6 _1 fc, ui — > b~ z ui, 
and w — > b^ +d+z w, where b = e l with I — > 0. The flow 
equations for the parameters T, A, e and v are obtained 
through a direct application of the RG transformations to 
eq. ([}]) in the Fourier space resulting, to one-loop order, 



dT 



cQ 2 ), 



de 

di 



2e, 



dv 
dl 



(z - Qv. 



(9) 



(10) 
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where c = j^P^T 1 ^- £ ' ^ d * s ^ e complete solid angle 
in d dimensions, and Q2 — J q A(q)q 2 is the second mo- 
ment of the noise correlator. Here we have implicitly 
assumed that TA 2 ^ e and TA 2 3> Xv/a±, as it is ex- 
pected near the critical point where both e and v should 
vanishes. The long-range nature of demagnetization ef- 
fects is, within our RG approach, irrelevant. The same 
result will be obtained if one considers a local demag- 



netization term — eif. Hence, all the analysis developed 
below will also hold for local demagnetization effects. 

The renormalization of the moments of the noise cor- 
relator (Q n — f A(q)q n ) is obtained considering a ver- 
tex function |§,|f . The renormalization equations for the 
moments can be put all together in a renormalization 
equation for the noise correlator A(if), obtaining 



dA(H) 
dl 



{£-2()A{H)+(HA , (H)-c 



dH 2 



-A 2 (if) 



A(2T)A(0) 



(11) 



We look for a fixed point solution A* (if) to this equa- 
tion, obtained setting dA/dl to zero and adjusting £. 
Integrating from —00 to 00 one obtains 



c 



£ 

3' 



(12) 



provided dHA*(H) 0. On the other hand, for 
if -> we obtain A* (if) w Q + Q X \H\ - \Q 2 H 2 with 



CQ2 



e-C 



(13) 



Then, imposing scale invariance in eq. (|10|) and using 
eqs. ( p^|) and ( |l3| ) it results that 



(14) 



The exponents £ and z thus result identical to those ob- 
tained for the constant force case J8|,[)| . 

On the other hand, eq. ( p^| ) implies that there is a 
characteristic lenght £ ~ t~ v beyond which scale invari- 
ance is lost due to demagnatization effects, and a char- 
acteristic velocity v c ~ e e (or field rate h c — ev c ) beyond 
which the system is supercritical, where 



(15) 



The line v c ~ e e devides the phase diagram (v, e) into 
two regions. For v > v c the noise is annealed and, there- 
fore, the model is supercritical as discussed above. On 
the contrary, for v < v c the model is subcritical for all 
e > 0. The critical state is obtained in the double limit 
e — » and v = h/e — ► 0. The scaling exponents 7 = 1 
and v = 1/2 are exact while £ and z are given by eqs. 
( |l2| ) and (0), respectively, up to the first order in e. 

In the subcritical regime the dynamics takes place in 
the form of avalanches, characterized by the avalanche 
size P(s) = s~ T f(s/s c ) and duration P(T) = 
T~ a g(T/T c ) distributions, where s c and T c are the 
avalanche size and duration cutoffs. In the subcritical 
state s c ~ e -1 /" '. At criticality s c ~ L D and £ ~ L, where 



D is the avalanche dimension, leading to the scaling rela- 
tion a = 1/Dv. Other scaling relations are obtained tak- 
ing into account that x — (s) and J dsP(s) = J dTP(T), 
which lead to 7 = (2 — t)/<t and (r — 1)D = (a — l)z, 
respectively. Taking 7, v and D as independent expo- 
nents we obtain the scaling relations for the avalanche 
exponents 



7 

Dv' 



1 



D 



71/ 



(16) 



On the other hand, for d < d c , the avalanche dimension 
and the roughness exponent arc related via D = d + £. 
Hence, we can compute r and a using the RG estimates 
of Q and z and the exact values of 7 and v. Above the up- 
per critical dimension d uc = 4, one obtain the mean-field 
exponents D — d uc = 4, r = 3/2 and a = 2. 

Before compare our exponent predictions with experi- 
ments and simulations, let us make some remarks about 
experimental measurements. The Barkhausen signal 
V(t) is the voltage produced from a pickup coil around a 
ferromagnet subjected to a slowly varying applied field. 
In the low-frequency limit the time scale for domain wall 
motion is much smaller than the time between jumps and, 
therefore, one may guarantee that each induced voltage 
jump corresponds with a single avalanche in the domain 
wall motion. A resolution voltage level Vr is defined, 
such that one can not resolve details below Vr. An el- 
ementary Barkhausen jump can thus be defined as the 
portion of the V(t) signal delimited by two subsequent 
intersections of the signal with the Vr line. With this 
definition, the duration T is simply the time interval be- 
tween these two subsequent intersections and the size s 
is the area delimited by V(t) and Vr between the same 
points. 

Evaluating the scaling exponents for d — 2, using the 
scaling relations in eqs. (|l6|), the RG estimates of £ and z 
and the exact values of 7 and v, we obtain r = 5/4 — 1.25 
and a — 10/7 w 1.43. These values can be compared 
with Barkhausen measurements in magnetostrictive ma- 
terials [|]J^J^1. Earlier measurements by Urbach et al 
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|ij gives r = 1.33 ± 0.10. More recently, Durin and 
Zapperi (DZ) [^| reported the more accurate exponents 
t = 1.28 ± 0.02 and a = 1.5 ± 0.1, which are close to 
our estimates. In the last case the experiments were per- 
formed under an applied stress. Positive magnetostric- 
tive materials show an increase in internal domain wall 
lengths under applied stress, which make possible the 
analysis of finite size effects in experimental results. Ac- 
cording to DZ measurements, the signal amplitude cut- 
off increases with the stress while the avalanche dura- 
tion cutoff decreases with stress, in such a way that the 
avalanche size cutoff s c is nearly independent on stress. 
The stress dependence of the signal amplitude cutoff was 
also observed in j| . We still do not understand the scal- 
ing of the distribution cutoffs with applied stress. How- 
ever, the absence of stress dependence in s c suggests that 
the system is in a subcritical state, where the correla- 
tion length is smaller than system size. Unfortunately, 
the demagnetization factor is kept constant in these ex- 
periments and, therefore, one can not determine if it is 
actually a control parameter. 

On the other hand, numerical simulations of the cel- 
lular automaton version of eq. ([!]) have also been per- 
formed The more accurate numerical estimates, re- 
ported by Durin and Zapperi are r = 1.26 ±0.04 and 
a = 1.40 ± 0.05, which are in very good agreement with 
our estimates. However, as in experiments, the demag- 
netization factor is kept constant. Numerical simulations 
of the same model, but including long-range dipolar in- 
teractions, clearly shows that the demagnetization factor 
is a control parameter We expect the same result 
when long-range dipolar interactions are absent. 

Hence, we strongly suggest to change the demagne- 
tization factor in both, experiments and simulations. In 
numerical simulations it is reduced to change a model pa- 
rameter, while it can be done in experiments changing, 
for instance, sample geometry. Changing the demagneti- 
zation factor will be a very good test to our model predic- 
tions. After that, one may analyze if the avalanche dis- 
tribution cutoffs are a consequence of finite size effects 
or just an evidence that the system is in a subcritical 
state. Based on numerical simulations, Bahiana et al Q] 
rule out the second possibility. However, they may have 
taken the demagnetization factor so small that the corre- 



lation length is smaller than system size. Decreasing the 
demagnetization factor the system will, according to our 
predictions, evolve to the subcritical state. 

We conclude that the long-range demagnetization field 
does not change the universality class from that of local 
demagnetization field. This result contradicts previous 
conclusion by Urbach et al , based on numerical simu- 
lations in one and two dimensions and Barkhausen noise 
measurements. According to them, the addition of a 
long-range demagnetization field results in self-organized 
criticality, while a local demagnetization field does not. 
However, in our opinion, their analysis only reveals that 
long-range demagnetization effects are more appropriate 
to describe the domain wall motion in magnetic materi- 
als. As in other numerical simulations and experiments, 
the demagnetization factor is kept constant. 
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